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Abstract. We study the Dirac equation with slowly varying external po- 
tentials. Using matrix-valued Wigner functions we prove that the electron 
follows with high precision the classical orbit and that the spin precesses 
according to the BMT equation with gyromagnetic ratio g = 2. 



1 Introduction 



The high precision experiments of the anomalous magnetic moment of the 
electron or the muon use in their interpretation the BMT equation for spin 
precession [0. The derivation of Bargmann, Michel, and Telegdi does 
not even mention quantum mechanics, however. Rather their arguments are 
based on relativistic covariance together with a classical notion of an intrin- 
sic angular momentum. Thus one might wonder how exactly the classical 
BMT equation is connected with the quantum mechanical spin precession. 
Properly speaking, one should start from QED and then deduce a BMT- 
like equation for the spin motion. Such a program looks rather difficult at 
present. But even if we take the one-particle relativistic Dirac equation as 
starting point, as we do here, the problem seems little explored and is thus 
still of interest. 

Let us first recall the BMT equation. One considers a charge, mass m, 
charge e, moving in prescribed, time-independent external electromagnetic 
fields E = — V0, B = V X A and governed by the Lorentz force equation 



di 



m'yvt = e{E{qt) + c ^Vt x B{q^)) . 



Here is the position and Vt 
7 = l/^/U^JvJcj^. According to BMT, for g 
s along the given orbit is determined by 



q^ the velocity of the particle at time t with 
2, the precession of the spin 



d 

dt' 



s X 



mc 



B 



,7 



1 + 7 



-c X E 



[1.2) 



For g ^ 2 the pref actors of B, v x E change and also a term proportional 
to v{v ■ B) appears. We refer to the textbook by Jackson, Section 11.11, for 
details 0. 

The same physical situation can be described by the relativistic one- 
particle Dirac equation in the semiclassical limit. This means that the po- 
tentials 0, A are slowly varying on the scale set by the typical extension 
of the electron/positron wave function. On a formal level it is convenient 
to introduce the dimensionless scale parameter e, e -C 1, and to take the 
potentials to be of the form 



(ex) , Aiex) 



;i.3) 
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The state of the electron/positron is now described by the four component 
spinor ip which evolves according to the Dirac equation as 

d 

ih—ip = Htp (1.4) 

with the Dirac hamiltonian 

H = c7o(7 ■ {-ihVx A{ex)) + mc) + e(j){ex) (1.5) 

^ c 

acting as a hnear operator on ^^=1 L'^iM^)- Here 70, 7 are the Dirac gamma 
matrices, where we follow the conventions in For the semiclassical limit 
one starts with a well localized wave function ip, having zero positron com- 
ponent, and evolves it unitarily according to 

= e-^^*/'^7/> (1.6) 

over such a long time that the response to the external forces is visible. With 
our choice ( |1.3| ) this means times of the order e~^. If the standard picture 
is correct, then over that time-span tpt should have a negligible positron 
component and should be well concentrated on the classical path determined 
by (|1.1|) . In addition, the average polarization should satisfy ( [L.2| ). This is 
exactly what we are going to show. 

Before entering in the details of our argument, it might be useful to recall 
the semiclassical limit of the Schrodinger equation 

^^|^=(-|^A + n^=^))^ (1-7) 

with the slowly varying electrostatic potential V. If the initial ip is well local- 
ized at the classical phase point {q, p), then ipt is well localized at {qt,Pt) with 
q = -^p, p = —W{q). This is an extremely well-studied chapter of math- 
ematical physics with higher order corrections and sharp error bounds 
available ||^. The Schrodinger equation is scalar. Put it differently, if we 
set the potential to zero, then the energy function of (1.7) is E{p) = i^p^ 
which is a scalar function. In somewhat more complicated cases, like an 
additional short scale periodic potential 0, internal degrees of freedom, or 
several particles, for zero external potential the energy will be matrix- valued 
as a function of a suitable quasimomentum. If this energy matrix has nonde- 
generate eigenvalues (= energy band functions), then in essence one is back 
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to the scalar case. Interference between bands is destroyed through rapid 
oscillations and the evolution in each band separately is semiclassical 0. 
In general, energy bands may cross, i.e. the eigenvalues may be degenerate 
along lower dimensional submanifolds. The wave function can then propa- 
gate into neighboring bands giving rise to intricate interference patterns § . 
In the case of the Dirac equation, the energy matrix is a 4 x 4 matrix and it 
has the two energy bands E±{p) = ±ca/ m^c^ + both two-fold degener- 
ate, which is at the origin for a scenario very distinct from the scalar case. In 
fact, in the semiclassical limit, the position and momentum become classical 
variables, whereas the spin remains fully quantum mechanical. 

It should be noted that from a purely abstract point of view the Dirac 
equation is the simplest example for a matrix-valued energy function. One 
would like to have at least two bands each of which should be degenerate. 
This makes four to be the minimal matrix dimension. 

We close the introduction with some brief comments on previous work. 
The most elementary approach is to consider average values of spin and posi- 
tion and to use what is known as Ehrenfest's theorem. For the Dirac equation 
the details have been carried through by Fradkin and Good 0. The more 
sophisticated WKB approach to the Dirac equation was initiated by Pauli 
llO| and completed by Rubinow and Keller The full semiclassical, van- 



Vleck-type propagator has been worked out recently by Bolte and Keppeler 



12| including a trace formula for the eigenvalues. The WKB method has the 



drawback of being local in time only. In the case of the Schrodinger equa- 
tion one understands how to continue beyond the caustics by adapting new 
coordinates in the so-called Lagrangian manifold P, |13]. For matrix- valued 
energy functions such a program has not been attempted yet. Instead, we 
will use here matrix-valued Wigner functions. The abstract theory for this 
approach has been developed in and applications to beam physics are 



discussed in 15 



2 Semiclassical limit of the Dirac equation 

In the Appendix we study abstractly the semiclassical limit for Schrodinger 
type equations with matrix-valued symbols. In this section we merely tran- 
scribe these general results to the Dirac equation. The emergence of the 
BMT equation will be explained in Section ^. 

In the following p denotes the classical momentum and we use q, p as a 
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pair of canonical coordinates. {■, ■} is the Poisson bracket. We define it also 
for matrix- valued functions A, B through 

{A,B} = VyA-VqB -VqA-VpB (2.1) 

which is again a matrix-valued function. Here ■ refers to the scalar product 
between the gradients. Note that the order of factors must be respected. In 
particular, 



{A^}^0 (2.2) 



in general. 

The Dirac matrix reads 



Hd = C7o(7 ■ (P - -Mq)) + mc) + e<P{q) , (2.3) 



compare with ( |1.5| ), which is regarded as a 4 x 4 matrix- valued function on 
the classical phase space, {q,p) G x M^. Hd is easy to diagonalize. The 
eigenvalues are 

h± = ±cpo + ecf) (2.4) 
and the corresponding two-dimensional eigenprojections 

P± = i (1 ± ^70(7 ■ (P - -A) + mc)) , (2.5) 

where po = \/vn?(? + (p — f ^)^- Thus 

HD = h+P+ + h_P_. (2.6) 

P+ is the electron and P_ the positron subspace. In addition, we need the 
polarization matrix along the orientation a,\a\ = 1, which is given by 

a ■ 5 = (7 ■ 5 + -a ■ (p - - A))75 (2.7) 



with 



g = a- \ (a • {p - -A))(p - -A) , (2.8) 

Po(Po + mc) c c 
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cf. 0. a ■ S has eigenvalues ±1, (a ■ Sy = 1, and commutes with Hd, 

[Hd,S]=0. (2.9) 

We discuss only the motion of an electron, the positron subspace being 
simply its mirror image. In the semiclassical limit the state of the electron 
in general is described through the matrix-valued Wigner function Wt{q,p) 
on classical phase space. According to the Appendix, Wt is governed by 

^Wt{q,p) = -{h+{q,p),P+Wt{q,p)} -t[Hi^\q,p),Wt{q,p)] (2.10) 
with 

//W = -z([P+,{/i+,P+}] 

+ ^P+(/l+{P+, P+} - /l_{P_ , P-})P+) = Hte + Hrrn ■ (2.11) 

If the initial wave packet is well localized at q^,p^, then one has to solve 
( |2.1CI| ) with the initial condition Wo{q,p) = 5{q — q^)6{p — p'^)\(po) {(po\ and 
(fo £ P+C^. In the course of time, this structure is preserved, 

Wt{q,p) = 6{q - q_MP " P-t)\Vt) {y^t\ , (2.12) 

where qt,Pt is the solution to 

q = {h+,q} = Vph+, p= {h+,p} = -Vqh+ (2.13) 

with initial conditions q^,p^- ( p.l3| ) is equivalent to ( |1 . 1| ) with po = mc-f and 
the velocity defined through 

v = —{cp-eA). (2.14) 

Po 

The spinor ipt evolves with a time-dependent hamiltonian according to 

i^^t = H+{qt,Pt)ipf (2.15) 

By ( pID this implies <^t e P+(qi,pJC^ for all t. 

In (|2.10|) the Poisson bracket ■} transports the spinor along the clas- 
sical orbit maintaining its orientation. Therefore, in general, the spinor would 
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move out of the electron subspace which is precisely compensated for by Hhg. 
For that reason, the analogue of Hf,f, has been baptized Berry term in |jl6 



while the remainder, if„n, has "no name" and is referred to also as Poissonian 



curvature in |jT7|. 



It remains to compute iji^"*. Hf,e is actually not needed in explicit form, 
but we list it for completeness, 

Hbe = 7^(75707 ■ {c-'v X F) - t-^f ■ F) (2.16) 

with F = E + c^^v X B. For the no name term, we note that 

e 

2~p'o 

and hence 



{P+, P+} = {1 - P_, 1 - P_} = -^^75707 ■ B (2.17) 



Hnn = ^«P+((-cpo-e0){P+,P+} + (-cpo + e0){P_,P_})n 

= --^P+757o7-BP+. (2.18) 

3 Spin precession 

In the semiclassical limit the average electron polarization along a is given 
by 

a-st = I d'qd'ptT[Wt{q,p)P+iq,p)a-Siq,p)P+iq,p)] (3.1) 

with Wt the matrix- valued Wigner function. Note that [P+, S] = 0. Accord- 
ing to the Appendix, we have 



d d 
—a ■ St = — 

dt dt 



j d^qd^pti[WtP+a- S] (3.2) 
^ j d\ d^ptT[P4{Wt, Hd] - {Hd, Wt])P+a ■ S] 
J d'q d^ptT[P+{{Wt, + t[Wt, Hnn])P+a ■ S] 
j d'q rf3ptr[W^i({/i+, P+a • 5P+} + P+a ■ 5P+])] . 



Using P+{/i+, P+}P+ = we have 

^a-st= [ d''qd'ptT[WP4{h+, a-S}- 2-^[757o7 ■B,a S])P+] (3.3) 
at J 2po 

and only have to compute the term P+(. . . )P+, which is accomphshed using 
the identities 

ipo/mc)P+v ■ 775P+ = P+v ■ SP+ = cP+^5P+ . (3.4) 
For the Berry term we have 

P+{/i+, a ■ 5}P+ = P+a ■ 5P+ = P+((7 ■ ^ + c'^a ■ vh5)P+ 
= c-'{- ^1 c~'{a . v){v . V) + -^^{a ■ v))P^v ■ 5P+ 

_^-2^2_^a.v)P^-fvP^, (3.5) 

mc + Pq 

where we used ( p.4| ). The last term transforms as 



P+7 • vP+ = P+v ■ SP+ — — c'\v ■ v)P+v ■ SP+ (3.6) 

mc + po 

and therefore 



P+a-SP+ = — — — c-\-(a-v)P+v ■ SP+ + (a-v)P+v ■ S P^ 

mc + Po 



Po 



c-^P+{a X (v X v)) ■ SP+ (3.7) 

mc + Po 

^^c-^P+a -{Sxivx v))P+ 

mc + Po 

^ ^ P+a ■ {S X {c-^v X E + c-^v X {v X B)))P+ . 



mc 1 + 7 

For the no name term we obtain, using again (^ 

-i^P+[^5lo-fB,a-S]P+ 

= -P+j ■ {B X g)^,P+ 
Po 

= -{P+{Bxg).SP+ ^^c-^v-{Bxg)P+vSP+) 

Po rnc + Po 

= —P+a-\-(SxB) + ^—c-^(Sx(vx(vxB)))]P+. (3.8) 

mc 7 1 + 7 
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We add Q and and insert in Eq. (^. Using (^U^ ) we finally obtain 

^a- s = —a- \s X (-B —c'^vxE)], (3.9) 

dt mc ^ ^-f 1 + 7 ^ ^ 

which is the BMT equation for g = 2. 



A Appendix 

We want to understand the solution of the Dirac equation with slowly varying 
potentials as in (|1.3|) . Changing to the rescaled space-time variables, denoted 
again by the same symbols, the Dirac equation becomes 

d 6 
ie/i— = [7o7 • {-iehV^ - -A{x)) + mc-fo + e(f){x)]ilj. (A.l) 

Since only the combination eh appears, we may set /i = 1 at the expense 
that formulas do not look right dimensionally. 

At this point, it is more convenient to abstract from the particular form 
of the Dirac equation and to consider a general unitary evolution governed 
by 

te-^^ij = H'ij (A.2) 

with tp an n-component spinor, ip G L"^ (M'^) ■ To define H"^ we start 

from a symmetric nxn matrix- valued function H{q, p) on the classical phase 
space r = M'^ X Mf^. Then, in spirit, H'^ = H{x, —ieVx)- While for the Dirac 
operator this substitution is unambiguous, in general it is not, and we define 
if^ through the Weyl quantization 

H^^ix) = (27^)-" I d'y I d'^H (^,^^) e^^--yy^i,{y) . (A.3) 

Next we have to list the properties of the symmetric matrix- valued Hamilto- 
nian function H. We assume that it has the spectral decomposition 

m 

H{q,p) = Y,h,{q,p)P,{q,p). (A.4) 
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hj are the real band energies. We order them as hj < hj^i and we assume 
that bands do not cross which means that 

hj{q,p) < hj+i{q,p) , (A.5) 

j = 1, . . . ,m — 1. If bands cross, then the situation becomes considerably 
more complicated. We refer to [|] for the case of two nondegenerate bands 
touching each other at a single point. Pj{q,p) are orthogonal, symmetric 
projections. Since H is assumed to be smooth in q,p and by ( |A.5| ), their 
degeneracy must be independent of q, p, 

tTPj{q,p) = dj (A.6) 

with degeneracy dj. We have Yl]jLi = n. dj = 1 is a scalar band. We study 
in the limit e —>■ through the evolution of Wigner density matrices. 



which for the pure state ipt is defined by 

W,%q,p) = {2n)-'l d'iMx - '-^WAx + \i)e^^-^ (A.7) 



and for a general density matrix by the corresponding incoherent superpo- 
sition. W^ is a n X n matrix-valued function on phase space. Using the 
evolution law (|A.2|) for ipt and expanding the right hand side of ( |A.7|) in e. 



one obtains 



= -e'H[H, Wl\ + \{{Wl H} - {H, W^}) + Rl (A.8) 



with a remainder Rl which vanishes as e ^ ||14|. Note that H is the given 
n X n matrix-valued hamiltonian function. In ( |A.8| ) the first piece of the 
generator is rapidly oscillating while the second piece is of order one. Such 
a situation has been studied abstractly by Davies [|T^. He proves that, for 
any t ^ 0, — Wt in the limit e — and that Wt commutes with H, 
[H, Wt] = 0. The order one piece of the generator is projected onto the 
invariant subspace of the evolution generated by —i[H, ■]. This means that 
in the limit e ^ ( [A.8| ) goes over to 

-Wt = J2P^^i{Wt,H} - {H,Wt})P, (A.9) 
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with initial conditions given through hnie^o SjLi Pj^o^j = ^o- Note that 
for all t we have [Pj, Wt] = 0, j = 1, . . . , m. It is instructive to rewrite (|A.9|) 
in a somewhat different form. We insert the spectral decomposition of H, 

^ m 

- J2 H} - {H, W})P, (A.IO) 

^ m m 

j=l i=i 

m -.mm 

= - E ^^^^ + 2 5^ 5Z h^mW, Pi} - {Pe, W})P, . 

j=l j=l 1=1 

To continue we need two identities. The first one reads 

PAh,,P,}Pj = (A.ii) 

which follows from {hj, Pj} = {hj, P^} = {hj, Pj}Pj + Pj{hj, Pj}. Therefore, 
using that PjW = WPj, 

P,{h,, W}P, = {h,, P,WP,} - P,W{h,, P,} - P,}WP, 

= {h,,P,WP,}-[P,WP,,[P„{h,,PM- (A.12) 
The second identity is derived from 

A{B, C} - {A, B}C = {AB, C} - {A, BC} , (A.13) 

compare with Then 

P,{P,,W}-{P,,P,}W = {P,,W}-{P,,P,W}, 
W{Pj,P,}-{W,P,}Pj = {WP,,Pj}-{W,P,} (A.14) 

and 

Pji{W, Pj} - {Pj, W})P, = WP,{P„ P,}P, - P,{P„ P,}P,W 

-PjiPAPj, W} - {Pj, W}Pj - {PjW, Pj} + {Pj, WPj})Pj 

= [W,P,{P„P,}P,], (A.15) 
where we use [W, Pj] = and again ( |A.13| ). In addition, again by (|A.13|) , 

Pl{P2, P2} = -{Pi, P2}{1 - P2) , {P2, P2}Pl = -(1 - P2){P2, Pi} ■ 

(A.16) 
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Then, by (^J^ and P1P2 = 



Pi{P2,W}-{P,,P2}W = -{P,,P2W}, 

W{P2,Pl}-{W,P2}Pl = {WP2,Pl} (A.17) 



and 



Pl{{W,P2} - {P2,W})Pi 

= PliW{P2, Pi} - {Pi, P2}W + {Pi, P2W} - {WP2, Pl}Pl 
= WP^il - P2){P2, Pl}Pl - Pl{Pl, P2}(1 - P2)PlW^ 

= -[iy,Pi{P2,P2}Pi], (A.18) 
where used that P^W = PiW, Pi(l - P2) = Pi, (|XT^) , and 

Pl{WP2,Pi} - {Pi, W^P2}Pi = {Pi,P2H^Pi} - {P,WP2,Pl} = 0, (A.19) 

which follows from ( [AA3|) , P2W = WP2, and P1P2 = 0. Let us define 

m 

Hi^^ = -t[P,,{h„P,}]-'-h,P,{P„P,}P,+ '- h,P,{P,,Pe}P,. 

(A.20) 

Then we have transformed ( |A.9| ) into 

g-Wt = J2{-{h,,P,WtP,} - t[Hi^\P,WtP,]). (A.21) 
i=i 

If we assume that the initial Wq is concentrated in band j, i.e. PjWoPj = 
Wq, and solve 

^Wt = -{h„Wt}-^[Hi^\W,], (A.22) 

then PjWtPj = Wt at all times. If Wq is a pure state for all q,p, Wo{q,p) = 
\i'o{<l,p)){4'o{Q,p)\ with PjV'o = "ipo, then this structure is maintained in the 
course of time, Pjipt = 'ipt, and ipt solves 

^^^^;, = -{h^,^P,}-^Hi^)^l;,. (A.23) 
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In particular, if initially Wo{q,p) = 6{q — q^)S{p — P°)|v^o)(v^o| with ipo e 
PjC", then at any other time 

Wtiq, p) = 5{q - q_^)5{p - p_,)\^t) , (A.24) 

where {qt,Pt) is the solution to 

q = {hj,q} = Vphj , p= {hj,p} = -Vghj (A.25) 

with initial conditions q^,p^. Thus the wave packet propagates along the 
classical orbit with the band function hj as hamiltonian. ipt evolves according 
to the time-dependent spin hamiltonian Hi''\t) = Hi^\q^,Pf) of the j-th 
band as 

Note that Pj(gj, pjv^t = 

The final answer ( |A.24| ), ( |A.26| ) has one rather surprising feature. We 
focus on the j-th band. To compute the q, p evolution we only have to know 
the band function hj of the band under consideration. The first term in 
the spinor hamiltonian just ensures that, as the wave packet propagates, the 
spinor does not move out of the j-th band subspace. Clearly, [P^-, {hj^Pj}] 
knows only about hj and Pj, as does the second term of Hs'^\ However, the 
third term depends on all the other bands. Thus if we modify H"^ keeping 
the j-th band fixed, then the spinor motion, and only it, will be modified. 
Physically, one might have been tempted to argue that the other bands are 
separated from the band j by a large energy and therefore modifying them 
should leave the dynamics in the j-th band unaffected. Our computation 
shows that the actual behavior is otherwise. 
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